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Abstract 

We consider the exact solution of a model of correlated electrons based on the super- 
algebra Osp{2\2). The corresponding Bethe ansatz equations have an interesting form. 
We derive an expression for the ground state energy at half filling. We also present the 
eigenvalue of the transfer matrix commuting with the Hamiltonian. 
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Recently supersymmetric generalizations of the Hubbard model have attracted consider- 
able interest for their possible relevance in describing correlated electron systems. A successful 
example is the Essler, Korepin and Schoutens extended Hubbard model due to its super- 
conductive properties. This model is based on the S'Z(2|2) superalgebra and is exactly solvable 
in one dimension by the Bethe ansatz [|l|, 0]. Other interesting system is a supersymmetric 
free-parameter model based on the continuous representation of the gl{2\l) superalgebra dis- 
cussed by Bracken et al 0. We recall that, for a special value of the coupling constant, this 
model was also considered by Karnaukhov |^ . Shortly afterwords, an anisotropic deformation 
of this model was found ^ as well as its solution by the coordinate Bethe ansatz approach 

The main feature of these models is that in one dimension they are integrable and therefore 
they provide us with non-perturbative information concerning physical properties. In fact, they 
can be derived from supersymmetric solutions of the Yang-Baxter equation invariant either by 
the S'/(2|2) or by the Osp{2\2) superalgebras [0, ||, ^ [1^. However, at least from the Bethe 
ansatz point of view, there exists a supersymmetric solution which has not yet been discussed 
in the literature so far. This solution was found sometime ago by Deguchi et al |]12| in the 
context of a g-deformed Osp{2\2) superalgebra and very recently has been interpreted as being 
invariant by the twisted Ug[Sl{2\2)^'^'^] superalgebra ||T3|. It seems that this model exhausts 
the cases which can be derived by exploiting the Osp{2\2) superalgebra [1^. This then 
provides us with an extra motivation for studying this supersymmetric system. The purpose 
of this paper is to discuss this solution in the context of an exactly solvable model of correlated 
electrons by presenting its Bethe ansatz solution. 

The one dimensional lattice Hamiltonian is derived from the logarithmic derivative of the 
i?-matrix solving the supersymmetric Yang-Baxter equation. For explicit expressions we refer 
to refs. [|12], ^ |13|. Its turns out that the corresponding Hamiltonian can be written in terms 
of fermionic operators c| ^ and q^o- acting on the sites i of a chain of length L and earring spin 
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index o" = ±. We found that the Hamiltonian reads as 

L 
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where nj.o- = cj^Ci^o- is the number operator for electrons with spin a on site i and we write = 
rii^^ This Hamiltonian conserves the number and of electrons with spin up and 

down, respectively. Later on we will use the total number of electrons A^e = ^+ + and the 
number as the good "quantum" numbers for charactering the spectrum of the Hamiltonian. 
We observe that in the derivation of expression (1) we have made use of convenient canonical 
transformations and also we have disconsidered terms ( such as the one proportional to nj — rij+i 
) which are automatically canceled out when periodic boundary conditions are imposed. 

The Hamiltonian (1) presents extra fine tuned hopping terms, on site and off site Coulomb 
interactions. We recall that such interactions resemble much those appearing in the model of 
hole superconductivity proposed by Hirsch ||I4|, and they can be derived from first principles. 



This model becomes integrable when the couplings Vi and V2 are constrained on the unitary 
circle, i.e. + ¥2=!- The parametrization of these couplings in terms of the g-deformed 
parameter q = exp(i7) (0 < 7 < tt) of the Osp{2\2) superalgebra can be read as 

V^i =sin(7),V^2 = cos(7) (2) 
From equation (1) we see that when 7 — > 0, this Hamiltonian reduces to the model proposed 



by Essler, Korepin and Shoutens This has been earlier noted in ref. ||T3[. It is still possible 
to reduce the interval for the anisotropy 7 if we perform particle-hole transformations. Indeed, 
the combined set of canonical transformations (I) c|,^ ^ Ci^a- and (II) c|,^ (— l)Vj,o- leads to 
conclude that the spectrum at certain values of 7 and 71 — 7 are related to each other. More 
precisely, we have the following identity 

i7(7) = -//(vr-7) (3) 



and consequently we can restrict our analysis only for the regime < 7 < 7r/2. This identity 
indicates that the symmetric point 7 = 7r/2 is somewhat special. In fact at 7 = 7r/2 we are left 
with "almost" a free-fermion theory. In this case the spectrum splits in two sectors, depending 
on the parity of the total number of electrons in the lattice. We have two decoupled XY 
models, with antiperiodic or periodic boundary conditions for N^. even or odd, respectively. 
We note that H{'k/2) corresponds also to the limit f/ ^ 00 in the Bariev chain [1^, which is a 



special case of a generalized XY model proposed long ago by Suzuki [|1^ . This means that the 
Deguchi et al's _R-matrix provides also embedding (" covering" vertex model ) for this peculiar 
limit of the Bariev chain . 



We now turn to the diagonalization of Hamiltonian (1) by the coordinate Bethe ansatz 
formalism. In other words, we would like to solve the eigenvalue problem, H^> = E{L)^> , 
provided the wave function in the sector of A^^e electrons distributed on the positions 1 < 
a^Qi < xq^ < • ■ ■ < xqj^^ < L has the following form 

(4) 

p j=i 

where the P summation extends over all the permutations of the momenta (Pi, ■ ■ ■ , PnJ aiid 
sgn is the sign of the permutation. When electrons are enough far apart {\xq- — xq.\ > 2) it 
is direct to derive that the eigenvalues are given by 

E{L) = 2 cos(pj) + V2{L- 2N,) (5) 
i=i 

The next step is to consider the matching condition for the Bethe ansatz wave function. 
Since the Hamiltonian (1) has been derived from a factorizable i?- matrix, we can restrict 
ourselves to the discussion of the two-body problem (two electrons in the chain). This is 
given in terms of the two-body S-matrix, which connects the scattering amplitudes between 
states {(pi, cTi); (p2, (^2)} and {(p2, (^'2)'-, (Pi? "^i)}- We found that the non-null two-body S'-matrix 
elements are given by 

SXX{\) = SZZ{\) = 1, SXZ{\) = 5I+(A) = . f'^^^^y S-t{\) = Sl-{\) - ^i^^(227) 



sinh(A + 2i7)' ' ' sinh(A + 2^7) 

(6) 



where A = Ai — A2 and Xj are the "dressed" momenta rapidities which are related to the 
momenta pj by the following relation 

Now we have the basic ingredients to derive the Bethe ansatz equations. The next step 
is to face the problem of diagonalizing the spins degrees of freedom which are encoded in the 
^-matrix (6). However, since this S'-matrix is of 6- vertex type, this later problem can be solved 



by standard algebraic methods ||T8[. In the course of the solution we have to introduce new 
spin rapidities fij, j = 1, ■ ■ ■ , It turns out that the dressed momenta and the spin variables 
satisfy the following nested Bethe ansatz equations 



sinh(Aj/2 - i7/2) 



sinh(Aj - yUfe - ij) 



sinh(Aj/2 + 27/2) sinh(Aj — fik + ^7) 

TT sinh(/xj - Afc - i'j) ^ _ tt sinh(/Xj - /Xfc - 2^7) ^. ^ ^ . . . 
sinh(/ij - Afc + i'y) ^J^ sinh(;Uj - + 227) ' ' ' + 

and, in terms of the rapidities Xj, the eigenvalues E{L) are given by 

^ cos(7) - cosh(Ai) 

The Bethe ansatz equations (8,9) have the unusual peculiarity that the "dressed" momenta 
variable Xj enters in different ways in the momenta (as Aj/2) and in the bare phase shift 
(right hand side of (8) ). This is a special feature which distinguishes this model from the 
free-parameter supersymmetric correlated electron system 0, ^, |3- Furthermore, near half- 
filling we have determined that while Xj are real roots, the spin variables fij form strings of 
the following type 

/^i = /^i + ^| (11) 
By substituting this ansatz in equations (8,9), taking their logarithmic and afterwards 
performing the thermodynamic limit L 00, we are able to obtain the integral equations for 
the densities a{X) and p(/i) for the variables Xj and fij, respectively. They are given by 

a{X) = + rf^4(A - ^, 7)p(;,) (12) 

ZTT J -00 



r+oa ^ 

27rp(/i) = / dXip'^ifi- X,2-f)p{X) + dXilj^i^i - X,-i)(y{X) (13) 



oo 



where 7) = 2 arctan[cot(7) tanh(x)] and ip2{X) 7) = 2 arctan[tan(7) tanh(x)] and the sym- 
bol ' stands for the derivative il)'i2{'^l d^l) = '^^^'^ix^"'^^ • These coupled integral equations can 
be solved by standard Fourier transforms. Considering the thermodynamic limit of E{L)/L 
and by taking into account the expression of cr(A), we find that the ground state energy per 
particle at half-filling is given by the expression 

. • / \ /""^ , coshb(7r/2 - 7)1 sinhfco-fTT - 7)1 , , ,^ ^, 

eoo = -4 sm 7 / duj . , \ + cos 7 14 

JO cosh[co'7r/2J smh[co'7rJ 

We would like to conclude with the following remarks. The experience we have gained 
solving the Hamiltonian (1) by the Bethe ansatz is also helpful if one wants to determine the 
eigenvalues of the transfer matrix of the underlying classical statistical model. Exact infor- 
mation about these eigenvalues is of considerable interest, because the transfer matrix is the 
generator of the many conserved currents commuting with the Hamiltonian (1). Considering 
that the eigenvalues A(A, {Xj,fij}) of the transfer matrix are analytical functions of the vari- 
ables {Xj,fij}, it is not difficult to start with an ansatz fulfilling this requirement. This is a 
phenomenological approach which goes by the name of analytical Bethe ansatz. Taking into 
account the Bethe ansatz equations (8,9) (analyticity conditions), the unitary and the crossing 
properties of the transfer matrix , we can find the following expression for the eigenvalues 

Aa XX 1^ , ,^,^^^ sinh(V2- A/2-27/2) , ,,,,,, ^ ^ cosh(A/2 - ^2 - ^7/2) 
A(A,{A.,/.,}) = [a(A)] n ,i,h(A./2-A/2 + .7/2) +[^^^^^ g ,,,h(A/2 - A./2 + .7/2) 

_ ^ Iff sinh(A/2 - A./2 + ^7/2) ^ sinh(/i, - A + 2^7) 
^ ^ \l\ sinh(A/2 - A,/2 - ^7/2) f}^ sinh(/x, - A) 

^ cosh(A/2 - A,/2 - ^7/2) ^ sinh(A - /i, + 2^7) \ 
y,cosh(A/2-A,/2 + z7/2)/i, sinh(A - /x,) ^ 

where the functions a{X), 6(A) and c(A) govern the behaviour of the transfer matrix on the pure 

ferromagnetic pseudovaccum ( this is equivalent to the totally empty or full band of electrons). 

For instance, if we take into account the i?-matrix elements of ref. []T2| we find that 

sinh(^7-A/2) _ tanh(A/2) _ sinh(A/2) 

^ ^ sinh(^7 + A/2)' ^ ^ tanh(^7 + A/2)' ^ ^ sinh(z7 + A/2) ^ ^ 
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and from equations (15,16) it is clear that all the residues of A (A, {A^, fij}) vanish, provided the 
variables {Xj,fij} satisfy the Bethe ansatz equations. We believe that this expression can also 
be derived more rigorously by using our recent formulation of the quantum inverse scattering 



method for Hubbard-like models 19 



In summary, we have solved a supersymmetric model of correlated electrons by the Bethe 
ansatz approach. The ground state energy at half-filling as well as the corresponding transfer 
matrix eigenvalues have been also determined. Preliminary analysis indicates that the excita- 
tions in the model are gapless near half filling. Finally, from the Bethe ansatz point of view, the 
limit 7 — needs to be performed having extra care. Due to the peculiar toplogy of roots fij 
(see equation (11) ), the situation is analog to the Su{3) isotropic limit of the Izergin-Korepin 



model |^0[- We also remark that when both Xj and fij are real roots, one easily recover the 



t — J sector present in the 5/(212) model [jT|, 0. 
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